Orbital Evolution of Algol Binaries with a Circumbinary Disk 
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ABSTRACT 

It is generally thought that conservative mass transfer in Algol binaries causes their orbits to 
be wider, in which the less massive star overflows its Roche- lobe. The observed decrease in the 
orbital periods of some Algol binaries suggests orbital angular momentum loss during the binary 
evolution, and the magnetic braking mechanism is often invoked to explain the observed orbital 
shrinkage. Here we suggest an alternative explanation, assuming that a small fraction of the 
transferred mass forms a circumbinary disk, which extracts orbital angular momentum from the 
binary through tidal torques. We also perform numerical calculations of the evolution of Algol 
binaries with typical initial masses and orbital periods. The results indicate that, for reasonable 
input parameters, the circumbinary disk can significantly influence the orbital evolution, and 
cause the orbit to shrink on a sufficiently long timescale. Rapid mass transfer in Algol binaries 
with low mass ratios can also be accounted for in this scenario. 



Subject headings: stars: general - 
matter 

1. Introduction 



An Algol binary is a semidetached binary sys- 
tem consisting of (1) an early type, main sequence 
primary component which does not fill its Roche 
lobe, and (2) a lobe-filling, less massive star that 
is substantially above the main sequence. The less 
massive star is cooler, fainter and larger (Giuricin 
& Mardirossian 1981; Peters 2001). It is believed 
that star (2) is initially more massive, and evolves 
first to overflow its Roche lobe to transfer mass to 
star (1). After the rapid mass exchange, the lobe- 
filling star (2) becomes less massive in the Algol 
binary stage. 

The evolution of Algol binaries has been inves- 
tigated extensively (e.g., Plavec 1968; Paczynski 
1971; Refsdal & Weigert 1969). Although conser- 
vative mass transfer seems to roughly reproduce 
the observed characteristics of a considerable frac- 



binaries: close — stars: mass loss — stars: evolution — circumstellar 



tion of Algol binaries (e.g. Nelson & Eggleton 
2001; De Loore & van Rensbergen 2004), noncon- 
servative evolution with mass and orbital angu- 
lar momentum loss taken into account is required 
for better comparison between theories and ob- 
servations (Thomas 1977; Refsdal et al. 1974; 
Sarna 1993). With regard to the orbital evolu- 
tion, conservative mass transfer in Algol binaries 
leads to increase of the orbital periods because 
mass is transferred from the less massive star to 
the more massive component (Huang 1963). How- 
ever, there is strong evidence showing that the 
orbital periods of some Algol binaries are actu- 
ally decreasing (Qian 2000a,b, 2001a,b,c, 2002; 
Lloyd & Guilbault 2002). Although it is uncertain 
whether the measured period changes are long- 
term (secular) changes or transient fluctuations, 
they clearly demonstrate that these systems must 
undergo mass and angular momentum loss during 
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the evolution. Moreover, Qian (2002) find that the 
decreasing rates Porb scale with the orbital periods 

Porb roughly as -Porb OC Porb- 

This orbital decay may be explained by angular- 
momentum loss via magnetic stellar winds (Ver- 
bunt & Zwaan 1981). Sarna, Muslimov, & Yerli 
(1997) and Sarna, Yerh, & Muslimov (1998) pro- 
posed that dynamo action can occur in mass- 
losing stars in Algol-type binaries, and produce 
large-scale magnetic fields, which lead to mag- 
netic braking of the stars. However, there could 
exist some potential difficulties for the magnetic 
braking mechanism (see below). Here we sug- 
gest an alternative interpretation assuming that a 
small fraction of the transferred mass would form 
a circumbinary (CB) disk surrounding the binary 
system (see also van den Heuvel 1994). Previous 
works (Taam & Spruit 2001; Chen et al. 2006) 
have already indicated that the CB disk can effi- 
ciently remove the orbital angular momentum and 
accelerate the binary evolution in accreting white 
dwarf or neutron star binaries. We describe the 
input physics that is necessary for the evolution 
model in section 2. Numerically calculated results 
for the evolutionary sequences of Algol binaries 
in two cases are presented in section 3. We make 
brief discussion and conclude in section 4. 

2. Input Physics 

We consider an Algol binary consisting of a 
lobe- filling, less massive star (of mass M2), and 
an early type, more massive companion (of mass 
Ml). The total orbital angular momentum of the 
system in a circular orbit is 



J = a /U 



2w 

Porh 



(1) 



where a = (GMTPorb/47''^)^/^ is the binary sep- 
aration and fi = M1M2/MT is the reduced mass 
of the binary system, G, and Mt are the grav- 
itational constant and the total mass of system, 
respectively. We neglect the spin angular momen- 
tum of the components because it is smaller com- 
pared with the total orbital angular momentum of 
the system. Similar to Rappaport et al. (1983), in 
the calculations we consider two types of mecha- 
nisms for orbital angular momentum loss from the 
binary system, which are described as follows. 



2.1. Magnetic Braking 

Vcrbunt & Zwaan (1981) point out that single 
low-mass main-sequence stars will undergo brak- 
ing by the coupling between the magnetic field and 
the stellar winds. The magnetic braking mecha- 
nism assumes that the ionized particles stream out 
along the magnetic field lines which come off the 
convection envelope of the stars. The loss of the 
specific angular momentum in the stellar winds is 
very large because the outflowing material is tied 
in the magnetic field lines to corotate with the 
stars out to a long distance (> 5 — 10 stellar radii) 
(Schatzman 1962). In an Algol binary, the angu- 
lar momcntiim loss of the evolved star by magnetic 
braking cause it to spin down. However, the tidal 
forces of the system will continuously act to spin 
the evolved star up back into corotation with the 
orbital rotation. The spin up takes place with the 
expense of the orbital angular momentum. Hence 
magnetic braking indirectly carries away the or- 
bital angular momentum of the binary system. 

The mass loss due to the magnetic winds during 
the mass transfer is generally very small, so we 
can assume that the total mass is conserved, and 
the orbital evolution is connected to the angular 
momentum loss rate and the mass transfer rate as 
follows 



orb 
^orb 



3J, 



MB 



J 



3M2 

Mo" 



(1-9), 



(2) 



where q = M2/M1 is the mass ratio, and Jmb is 
the angular momentum loss rate due to magnetic 
braking. 



Jmb 



-3.8 X lO-^°M2P0(i?2/P0)^a;^dyncm, 

(3) 

where 7 is a dimensionless parameter in the range 
from to 4, and uj is the angular velocity of the 
star (Rappaport et al. 1983). Since q < 1 and 
M2 < 0, the second term on the right hand side of 
Eq. (2) always causes the orbit to increase. Orbital 
shrinkage occurs only when the magnitude of the 
first term is larger than that of the second one, 
i.e., 

-Porb ^ -^JmB 
Porb J 

One can show that with the magnetic braking law 
of Rappaport et al. (1983), the orbital decay rates 
with Eq. (4) are several orders lower than the ob- 
served ones. Moreover, from Eqs. (3) and (4) we 
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can derive the following relation, 



-Porb oc P' 



(27-7)/3 



orb 



(5) 



where we have used the relation for the radius of 
the lobe-filling star R2 = 0A62[q/{l + q)]'^/^a with 
0-1 ~ 9 ^ 0-8 (Paczynski 1971). As we can see, 
only when 7 > 4 can we have a positive correla- 
tion between — -Poib and Porb- More recent stud- 
ies based on observations of rapidly rotating stars 
show that the above magnetic braking form may 
be inadequate (Queloz 1998; Sills et al. 2000; An- 
dronov, Pinsonneault, & Sills 2003). Sills et al. 
(2000) instead propose an empirical angular mo- 
mentum loss prescription 



Jmb = —K-u,u ( — 

1712 



1/2 



T 2 / ^2 

•Jmb — — -f^w'^t^crit 

m2 



1/2 



W > Wcrit (6) 



where K = 2.7 x 10 gem , o^crit is the critical 
angular velocity at which the angular momentum 
loss rate reaches a saturated state, r2 and m2 are 
i?2 and M2 in solar units, respectively. Kim & De- 
marque (1996) suggest that Wcrit is inversely pro- 
portional to the convective turbulent timescale in 
the star at 200 Myr age. 



Wcrit — '^crit,0" 



T0 



(7) 



The magnitude of Jmb under this magnetic brak- 
ing law is at least one order of magnitude smaller 
than that in the Rappaport et al. (1983) model 
with 7 = 4 (Andronov et al. 2003), suggesting 
that more efficient angular momentum loss mech- 
anism may be required to account for the orbit 
decay in Algol binaries. 

2.2. CB Disk 

Once Roche lobe overflow occurs in semide- 
tached binaries, some fraction of the transferred 
material may leave the system in various ways. 
The mass outflow in cataclysmic variables (CVs) 
was noticed by Warner (1987) from the studies of 
UV and optical resonance lines. 

Mass loss can also lead to orbital decrease in 
mass exchanging binaries. We first consider the 
case that a fraction S of the transferred matter 
is lost from the binary system, carrying away the 



specific angular momentum of the secondary star 
(Rappaport et al. 1983). The change rate of or- 
bital period can be written as 



P 
P 



M2 
5— 
Mo 



l-q{l-5) 



1 5q 



l + q 3 1-hg 



(8) 

To account for the secular decreasing of orbital 
period, the mass loss fraction 5 from the binary 
system should satisfy 



5> 



3(1 -g^) 
3 - 2g - 3g2 



(9) 



For q = 0.5, it requires 5 > 1.8. Apparently 
this is impossible. Hence angular momentum loss 

via simple mass loss cannot interpret the orbital 
shrinking of some Algol binaries with low mass ra- 
tio. 

In a theoretical view, Spruit & Cao (1994) pro- 
posed that part of the outflow may be in the form 
of a slow wind near the orbital plane. It has been 
argued that during mass exchange the lost mat- 
ter may form a disk structure surrounding the bi- 
nary system rather than leave the binary system 
(van den Hcuvcl & dc Loorc 1973; van den Hcuvel 
1994). More recently, Taam & Spruit (2001) and 
Spruit & Taam (2001) suggested that a Keplerian 
CB disk could be formed as a result of mass out- 
flow in CVs. The CB disk can extract the orbital 
angular momentum from th binary orbit through 
tidal torques. 

Since the mass transfer processes are similar in 
CVs as in Algol binaries, we assume that a con- 
stant fraction 5 of the transferred mass feeds into 
the CB disk surrounding the Algol system, and 
Ml = — (1 — 5)M2 (if the mass transfer rate is 
sub-Eddington) . Same as in the standard thin ac- 
cretion disks, the viscous torque in the CB disk is 
(Shakura & Sunyaev 1973) 



T = -2'Kryllr—r, 
ar 



(10) 



where v is the viscosity, r the radius, S the surface 
density, and n = [G(Mi -|- M2)/r^]^/^ the local 
orbital angular velocity in the disk. 

Assume that the gravitational interaction of the 
binary with the disk occurs locally at the inner 
edge Tj of the disk^ . angular momentum from the 



^In the following we use the subscript i to denote quantities 
evaluated at 
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binary feeds into the disk at a rate proportional to 
the surface density Si. The viscous torque exerted 
by the CB disk on the binary can be shown to be 
(Taam & Spruit 2001; Spruit & Taam 2001) 



-'orb / \ ^vi 



(11) 



where 7^ = rj/a, t is the mass transfer time, and p 
is a parameter determined by the viscosity. Sup- 
pose u is the function of r and S with the following 
form (Spruit & Taam 2001) 



r \ / S 



(12) 



thenp= (m-|-l)/[2(2m-|-2 — n)]. In the following 
calculations we set m = n = 1, so p = 1/3. 

By simple algebra, from the above equations we 
obtain the angular momentum loss rate via the CB 
disk 



1/3 



where 



.7CB = 7- 



J 
A*' 



(13) 



(14) 



is the specific orbital angular momentum of the 
disk material at n. 

The viscous timescale tyi at the inner edge of 
disk is 

4r? 

*. = ^, (15) 

and Vi is given with the standard a prescription 
(Shakura & Sunyaev 1973) 



fi = assCsHi, 



(16) 



where ags, Cg, Hi are the viscosity parameter (in 
the following calculations we set ass = 0.01), the 
sound speed, and the scale height of disk at n, 
respectively. Assuming that the disk is hydro- 
statically supported and geometrically thin with 
Hi/n ~ 0.03 (Belle et al. 2004), we have 



(17) 



If we consider angular momentum loss only via 
the CB disk, the changing rate of the orbital pe- 
riod is then 



^ = -3- 



M2 

Mo 



1 - 7r?(l + q)- q{l 
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+ q 
(18) 



where rj ^ S^t/tyiY^^. Since M2 < 0, the orbital 
period will decrease when 



1 - g(l - (5) - (l/3)(5g/(l + g) 
7(1 + 9) 



(19) 



For conservative evolution, S = r] = 0, Eq. (19) 
recovers to g > 1 as expected. 

In Fig. 1 we plot 77 against q in the solid curve 
when the orbital period is constant and 5 <C 1. 
Here we take n/a = 7^ = 1.7 (Artymowicz & 
Lubow 1994). Figure 1 shows that the CB disk 
presents a plausible mechanism removing the or- 
bital angular momentum to explain the decrease 
of the orbital periods for reasonable values of t] (or 
5). Generally a larger r] is required for a smaller 
q. The dotted and dashed lines show two exam- 
ples of Algol binary evolution, revealing how the 
orbital period evolution depends on the initial pa- 
rameters. As seen from Eq. (14), the larger orbit, 
the larger angular momentum loss rate via the CB 
disk. This naturally results in a positive correla- 
tion between _Porb and Poih- 

3. Numerical Results 

In this section we present the results of numer- 
ical calculations of the Algol binary evolutions in 
two typical cases. We have adopted an updated 
version of the evolution code developed by Eggle- 
ton (1971, 1972) (see also Han et al. 1994; Pols et 
al. 1995). We consider both the angular momen- 
tum loss mechanisms described in section 2, adopt- 
ing the modified magnetic braking law (Eqs. [6]- 
[7]) for evolved stars with mass < 1.5Mq. We take 
solar chemical compositions {X ~ 0.7, Y = 0.28, 
and Z = 0.02) for both stars, and the ratio of the 
mixing length to the pressure scale height to be 
2.0. 

(1) Initial masses (M2,Mi) = (3.O,2.O)M0 and 
orbital period Porb.i — 1-4 days 

In a proto- Algol system with initial orbital pe- 
riod Porb,i <~ 1 — 2 days, the more massive 
star will fill its Roche lobe when it is still on 
the main sequence. This kind of Case A binary 
evolution has been investigated extensively (e.g., 
Eggleton 2000). Here we calculate the evolution- 
ary sequences for a binary with initial masses 
M2 = 3Mq, Ml = 2Mq, and the initial orbital 
period Porb.i — 1-4 days, assuming that a small 
fraction 6 (= 0.01, 0.005) of the transferred mass 
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forms a CB disk surrounding the system. Figure 
2 shows the mass ratio, the orbital period and the 
mass transfer rate as a function of the time of mass 
transfer. It is clear that the mass transfer proceeds 
initially on short (< 1 Myr), thermal-timescale, 
then on a much longer (hundreds of Myr) nuclear- 
timescalc when q < 1, this is why we usually ob- 
serve less massive donor stars in Algol binaries. 
The orbital period first decreases until q ~ 0.6. Its 
further evolution depends on the effect of orbital 
angular momentum loss by the CB disk. When 
6 = 0.005 the evolution is similar as the conserva- 
tive evolution with continuously increasing orbital 
period. When 6 = 0.01 the angular momentum 
loss is strong enough that at q ~ 0.25, Porb turns 
to decrease. The second orbital decay phase oc- 
cupies most of the mass transfer lifetime. In both 
cases the mass transfer rates have a peak values 
^ 3 X lO^^M.Qyr^^ at the beginning of Roche lobe 
overflow, then last a plateau phase ~ 100 — 200 
Myr. It is also noted that, in both cases the mass 
transfer rates reach a peak when i ^ 150 — 300 yr, 
indicating that the CB disk-induced angular mo- 
mentum loss begins accelerating the mass transfer. 
(The evolution of tj was plotted against the mass 
ratio q with the dotted curve in Fig. 1.) 

(2) Initial masses (M2,Mi) = (3.0,2.0)Mq and 
orbital period Porb.i — 3.4 days 

The calculated results for a relative wide sys- 
tem with Porb.i - 3.4 days and S = 0.032, 0.015 
are presented in Fig. 3. Compared with those in 
Fig. 2, the mass transfer rates are much higher, 
and the evolutionary timescalc is smaller by a fac- 
tor more than 100, suggesting that these systems 
arc less likely to be observed. For both values of (5, 
the orbital period shows a transition from decrease 
to increase when q ^ 0.2 — 0.4. The phases of orbit 
shrinkage last about 20% of the total mass trans- 
fer lifetime. When S = 0.015, the mass transfer 
rate somewhat fluctuates during the late stages of 
the evolution though we have smoothed the mass 
transfer rate for clarity. 

The calculated evolutionary tracks in the HR 
diagram are plotted in Fig. 4. The bold curves 
denote the phases in which the orbital period de- 
creases. For comparison with observations, we also 
show the locations of four Algol systems which 
have decreasing orbital periods and well-measured 
parameters summarized in Table 1. As shown in 
the figure, the calculated results are roughly in 




Fig. 1. — The figure of rj against the mass ratio 
q. The solid curve denote the division between pe- 
riod increase and decrease. The dashed and dot- 
ted curves correspond to the numerical results of 
Algol binary evolution with the initial orbital peri- 
ods Porba = 3.4 d{S = 0.032) and 1.4 d (5 = 0.01) 
and initial masses (M2,Mi) = (3.0, 2.0)Mo, re- 
spectively. 
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Fig. 2. — The evolution of the mass ratio q (left), 
the orbital period Porb (center), and the mass 
transfer rate M (right) with the mass transfer time 
for a Algol binary with the initial orbital period 
Porb.i = 1.4 d. The solid and dashed curves cor- 
respond to 6 = 0.01 and 0.005, respectively. The 
mass transfer rate have been slightly smoothed for 
clarity. 




Fig. 3. — Same as Fig. 2, but for the initial or- 
bital period Porb.i = 3.4 d. The solid and dashed 

curves correspond to S = 0.032 and 0.015, respec- 
tively. The mass transfer rate have been slightly 
smoothed for clarity. 
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accordance with the observations, suggesting that 
they may represent the possible evolutionary se- 
quences that lead to the formation of these sys- 
tems. Figure 5 compares the calculated results of 
Algol binary evolution with the observational data 
in the Tes — Porb plane. It seems that our model is 
also compatible with the observed spectral types 
of the four systems. 

4. Discussion and Conclusions 

It is an open question why some of the Algol 
systems show orbital shrinkage while others not. 
Although it is generally believed that there must 
be orbital angular momentum loss during the evo- 
lution of these binaries, the appropriate mecha- 
nism(s) for angular momentum loss has not been 
verified. Based on theoretical model of Spruit & 
Taam (2001) and Taam & Spruit (2001), we sug- 
gest that a CB disk may play an important role in 
determining the orbital change in Algol systems. 
Note that our CB disk hypothesis does not ex- 
clude the effect of (magnetized) mass loss. On the 
contrary, the formation of the CB disk may be 
closely related to the mass loss processes in the 
binary evolution. The existence of the CB disks 
may be revealed by their infrared radiation, as in 
GG Tau (Roddier et al. 1996). The Monte-Carlo 
simulation by Wood et al. (1999) has reproduced 
the morphology of the near IR radiation which is 
in accord with the observations by Roddier et al. 
(1996). 

In this paper we have calculated the evolution- 
ary sequence for Algol binaries with typical masses 
and orbital periods, taking into account both mag- 
netic braking and the CB disk. The detailed calcu- 
lations indicate that the orbital evolution of Algol 
binaries can be significantly affected by the CB 
disk, especially for Case A mass transfer. With 
adequate values of 5, it is possible to account for 
the decrease of the orbital periods in some Algol 
binaries. In addition, the existence of CB disk 
can accelerate the evolution and enhance the mass 
transfer rates, as already noticed in the evolution- 
ary study of CVs (Taam & Spruit 2001). This 
suggests a plausible explanation for the rapid mass 
transfer observed in a few Algol binaries with very 
low mass ratios (Qian et al. 2002). 

Finally, we note that the magnitude of the 
mass feeding parameter 6 is highly uncertain, let 
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Fig. 4. - The calculated evolutionary tracks for 
Algol binaries in the HR diagram. The solid and 
dashed curves correspond to the initial orbital pe- 
riod Porb,i = 1-4 d (5 ^ 0.01) and Porb.i = 3.4 d 
{6 = 0.032), respectively. The bold curves denote 
the phases in which the orbital period decreases. 
The filled rectangles represent the locations of the 
four Algol binaries listed in Table 1. 
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Fig. 5. — Evolutionary sequences in the donor 
star's effective temperature Teg versus orbital pe- 
riod Porb diagram. The solid and dotted curves 
correspond to the initial orbital period Porb.i = 1-4 
d ((5 = 0.01) and Porb.i = 3.4 d{5 = 0.032), respec- 
tively. The bold curves denote the phases in which 
the orbital period decreases. The filled rectangles 
represent the locations of the four Algol binaries 
listed in Table 1. 
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alone its variation with the mass transfer rates. 
These uncertainties make it difficult to present 
direct comparison between observations and the- 
oretical predications for Algol binaries. How- 
ever, our analysis provides a reasonable mecha- 
nism for angular momentum loss in semi-detached 
binaries including Algol binaries. More detailed 
multi-wavelength observations can provide strin- 
gent tests for the CB disk model, and the theories 
of the evolution of Algol binaries. 

We thank the anonymous referee for his/her 
helpful comments improving this manuscript. 
This work was supported by Natural Science Foun- 
dation of China under grant number 10573010. 
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Table 1 

Observed Parameters of four Algol System 



Source 


P (d) 


P(10~'' dyr-^) 


Spectrum 


Ml (Mo) 


M2(Mo) 


logLi(L0) 


logL2(io) 


References 


X Tri 


0.9715397 


-0.142 


A2/G3 


2.3±0.7 


1.2±0.3 


1.21±0.06 


0.45±0.18 


1 


AT Peg 


1.146082 


-0.638 


A4/G 


2.22±0.065 


1.05±0.025 


1.19±0.025 


0.38±0.07 


2 


AF Gem 


1.2435012 


-0.0992 


B9/G0 


3.37±0.11 


1.155±0.038 


1.78±0.09 


0.75±0.03 


3 


TX UMa 


3.0632881 


-0.713 


B8/G0 


4.76±0.16 


1.18±0.04 


2.30±0.04 


1.17±0.065 


4 



"■References: 1 - Mezzetti et aL (1980), 2 - Maxted et al. (1994), 3 - Maxted & Hilditch (1995), 4 - Maxted et al. (1995). 



